Circuit quantum electrodynamics 
with a nonlinear resonator 



One of the most studied model systems in quantum optics is a two-level atom 
strongly coupled to a single mode of the electromagnetic field stored in a cavity, a re- 



search field named cavity quantum electrodynamics or CQED (Haroche and Raimond, 



2006). This extremely simple quantum system has nevertheless nontrivial quantum 
dynamics described by the Jaynes-Cummings Hamiltonian (Jaynes and CummingsJ 
[T963) . In its implementation at microwave frequencies, it has allowed the observation 
of many basic concepts of quantum mechanics such as the quantum jumps of the 
electromagnetic field (Guerlin et al., 2007), as well as the generation and tomography 
of nonclassical states (Deleglise et al., 20081. In the context of quantum informa- 
tion, elementary quantum gates have been realized using the cavity as a catalyst for 



entanglement between atoms passing successively through it ( Rauschenbeutel et al., 
[T999j). CQED has recently received renewed attention due to its implementation with 
superconducting artificial atoms and coplanar resonators in the so-called circuit quan- 
tum electrodynamics (cQED) architecture (Blais et al., 2004 Walbaff et al., 2004). In 
cQED, the couplings can be much stronger than in CQED due to the design flexibility 
of superconducting circuits and to the enhanced field confinement in one-dimensional 
cavities, compensating the shorter coherence times of superconducting qubits. This en- 



abled the realization of fundamental quantum physics (Hofheinz et al., 2009 Ansmann 



et al., 2009 



Palacios-Laloy et al., 2010) and quantum information processing (DiCarlo 



et al., 2009) experiments with a degree of control comparable to that obtained in 
CQED. Even though the physical implementation is different, the system is described 
by the same Hamiltonian in cQED as in CQED. 

The purpose of this chapter is to investigate the situation where the resonator 
to which the atom is coupled is made nonlinear with a Kerr-type nonlinearity, caus- 
ing its energy levels to be nonequidistant. The system is then described by a non- 
linear Jaynes-Cummings Hamiltonian. This considerably enriches the physics since 
a pumped nonlinear resonator displays bistability (Siddiqi et al., 2005), parametric 



amplification ( Castellanos-Beltran and Lehnert, 20071, and squeezing ( Castellanos- 



Beltran et al., 2008). The interplay of strong coupling and these nonlinear effects 



constitutes a novel model system for quantum optics that can be implemented exper- 
imentally with superconducting circuits. 

This chapter is organized as follows. In a first section we present the system con- 
sisting of a superconducting Kerr nonlinear resonator strongly coupled to a transmon 
qubit (Koch et al, 2007). In the second section, we describe the response of the sole 
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Fig. 1.1 Scheme of the nonlinear resonator. A Josephson junction of critical current Jo is 
embedded in the middle of a A/2 resonator. It is coupled to a 500 transmission line through 
a coupling capacitor Cc and probed in reflection by a microwave field. 

nonlinear resonator to an external drive, with a particular emphasis on the bistable 
regime. Our main interest is the probability that the resonator switches from its initial 
dynamical state to a state of higher amplitude under driving by a microwave pulse. 
We present here both measurements and numerical simulations of this process and 
obtain quantitative agreement. In the third section, building on our understanding of 
the switching process, we show how the resonator bistability can be used to perform 
a high-fidelity readout of the transmon qubit, a crucial task for quantum information 



processing (Mallet et ai, 2009). In the course of this readout, the nonlinear resonator 



is energized by a small microwave field; it is then natural to investigate what is the 



quantum backaction exerted by the intracavity field on the qubit (Ong et al., 2011). 
This is the subject of the fourth section. 

1.1 Presentation of the system 

1.1.1 Nonlinear resonator 



The heart of our device is a superconducting coplanar waveguide resonator ( Frunzio 



et al., 2005 Goppl et al., 2008) made nonlinear with a Josephson junction. Linear 
coplanar waveguide resonators consist of a section of length L of a coplanar waveguide 
transmission line, supporting modes at frequencies (D„ = (n + l)7rc/L, c being the 
speed of light in the transmission line. In the limit where the quality factor of each 
mode is high, one can consider only the dynamics of the mode closest in frequency 
to the system of interest (namely the qubit to which it is coupled), in our case the 
fundamental mode n = 0. The resonator Hamiltonian is then H = hwra^a, where 
Ldr — Qq and a (a^) is the usual annihilation (creation) operator. 

In our experiment a Josephson junction of critical current Iq is inserted in the 
middle of the transmission line resonator (see Fig. |1.1[ ). A Josephson junction can be 
considered from an electrical engineering point of view as a lumped dissipationless 
nonlinear inductance. As a result, the mode structure of the resonator is deeply mod- 
ified. First, the frequencies of all the resonator modes are differently shifted by the 
introduction of the lumped inductance, so that they don't follow the simple relation 
a)„ — {n + l)nc/L anymore. Second, the nonlinearity of the Josephson inductance 
results in nonlinear dynamics for each resonator mode and couplings of the resonator 
modes with each other. This could make the dynamics of this distributed nonlinear sys- 
tem hardly tractable. However, in the limit where the frequencies of the signals driving 



the resonator do not generate harmonics that are resonant with higher modes ( Zakka- 



Bajjani et al., 2011 ), one can approximate the resonator to only one nonlinear mode. 
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This is valid provided the frequency of higher modes is sufficiently different from 
{n + 1)1:0/ L (which requires a relatively large lumped element inductance) and pro- 
vided the frequency of the signals stays close to u)r (which requires that the quality 
factor of the fundamental mode is sufficiently large). All the results presented in the 
following are obtained within this approximation. More details on the modelling of the 



distributed nonlinear resonator can be found in (Wallquist et ai, 2006 



We now derive the Hamiltonian of the nonlinear resonator restricted to the funda- 



mental mode, and find explicit expressions for the nonlinear constants ( Palacios-Laloy 



et ai, 2008 Ong et al., 2011). The first step is to map the distributed nonlinear 
resonator onto an equivalent series combination of a lumped element inductance L^, 
capacitance Cc (which includes the junction capacitance) and Josephson junction of 



critical current Iq, as shown in Fig. 1.2 ^Manucharyan et al. , 2007 1 . In the limit where 



the Josephson inductance Lj is completely negligible compared to the resonator in- 
ductance, one can show that Lc = itZo/2uji and Ce = 2/7r.ZoWi, where Zq 50 fl 
is the resonator characteristic impedance and wi the resonance frequency in absence 
of the junction. When Lj is not negligible, one has to adjust numerically and Cg 
so that the impedance of the equivalent circuit fits the impedance of the distributed 
resonator. This is the approach that will be used in the following. 



Using the notation introduced in Fig. 1.2 we obtain the equivalent circuit Hamil- 
tonian 

where ipo = h/2e is the reduced superconducting flux quantum. Since the current / 
flowing through the inductance is the same as that in the junction, we also have 

/= — =/osm , (1.2) 

yielding an implicit relation (pi = g{4>) between the two phases. Eliminating (pi and 
expanding g{(j)) in powers of 0, we obtain the nonlinear resonator Hamiltonian to any 
order of Josephson junction nonlinearity. For instance, to sixth order we find 

2Lt 2Cc 24 LtfQ 720 LtfQ 

where Lt = £j + Lc is the total inductance and p = Lj / Lt the participation ratio of 
the Josephson inductance in the total resonator inductance. 

This Hamiltonian can be writen in terms of and a with <j) — iyJ~hZ^j2{a — a)) and 
q = yJh/2Z^{a + a' ) with = \/ L^_/Cc- Once expanded, the nonlinear terms 4>^ and 
(p^ yield products of creation and annihilation operators to various powers. Using the 
rotating-wave approximation, we keep only those with equal number of annihilation 
and creation operators yielding the Hamiltonian of a Kerr nonlinear resonator (KNR) 

K 2 K' 3 

Hnl/^^ u)ra)a^ a) ^ a) a^, (1.4) 

2 3 

where 
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Fig. 1.2 Equivalence between tiie distributed nonlinear resonator and a series combination 
of an equivalent inductance L^, capacitance Co and Josephson junction of critical current /q. 

K^-"^, (1.5a) 
2 K'^ 

K' = - (lOp-9). (1.5b) 

are the nonlinear constants respectively to third (Kerr constant) and fifth order, = 
Xj \jLtCc and Rk = h/e^. Here and below, we have dropped the resulting small 
correction to w^- Note that we will only consider ^ in the last section of this 
chapter where it should be taken into account for a quantitative agreement between 
theory and the data. 

In order to allow measurement of the resonator without perturbing its dynamics, 
the resonator is connected to a 50 O measurement line through a small coupling ca- 
pacitor Cc that determines its quality factor Q, or damping rate n = uJr/Q.H can be 
probed by a microwave signal of frequency ujp and amplitude Cp, which is reflected at 
the resonator input [see Figs. |l.l| and 1.6 a)]. The amplitude and phase of the reflected 



signal is measured by homodyne detection at room temperature after being amplified 
by a cryogenic amplifier. More generally, driving of the resonator is described by the 
Hamiltonian term 

Hd/h^ e<ie-'"''*a^+h.c., (1.6) 

d£{s,p} 

where two different drives have been considered: the pump drive (amplitude and 
frequency ujp) which will describe driving of the resonator close of its resonance fre- 
quency and the qubit drive (amplitude tg and frequency ujs) will describe driving of 
the qubit close to its resonance frequency through the resonator. 

The complete resonator dynamics is described in the usual Markov and rotating 
wave approximations, and for sufficiently weak nonlincarity K' <C Wr, by the master 



equation (Walls and Milburn, 2008 1 



p = CrP = [Hnl + Hd, p] + K(nth + l)2?[a]/0 + K.nt\j:>[a^]p, (1.7) 
with the dissipator 

V[A]p=]^{2ApA^ -A^Ap- pA^A), (1.8) 

and where rith = l/[exp{tiiOr/kBT) — 1] is the number of thermal photons in the bath 
coupled to the resonator at temperature T, and Ub is the Boltzmann constant. 
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Fig. 1.3 Circuit model of a split Cooper-pair box made of a superconducting loop including 
two Josephson junctions of Josephson energy Ejo threaded by a flux $, shunted by a capacitor 
C, and charge-biased by a voltage source Vg through a gate capacitor Cg yielding a charge 
bias iVg = CgVg/{2e). 

1.1.2 Qubit 

The artificial atom used in our experiment is a split Cooper-pair box operated in the 



phase regime (i.e. a transmon qubit, see (Koch et al, 2007 Schreier et ai, 2008 
Houck et al., 2009)). As illustrated in Fig. |1.3[ it consists of two superconducting 



islands connected by a superconducting loop (a SQUID) made of two Josephson junc- 
tions of Josephson energy Ejo and shunted by a capacitor C. The SQUID is threaded 
by a flux $ and charge-biased by a voltage source Vg through a gate capacitor Cg. Its 
Hamiltonian is H = Ec {N — Ng)^ — Ej{^) cos {9), N being the number of Cooper 
pairs transferred from one island to the other, Ng = CgVg/{2e) the charge bias, 
the gauge-invariant phase difference between the superconductors, Ec ~ (2e)^/2C 
the box charging energy and i?j($) = 2Ejq |cos(7r<I>/$o)| the SQUID flux-dependent 
effective Josephson energy. This Hamiltonian can be diagonalized yielding the qubit 
eigenstates {\i)} with frequencies {wj}. In this eigenbasis, the Cooper-pair box Hamil- 
tonian restricted to its M first eigenstates is 

M-l 

Hg/h= Y,^^^^,^ = '^u^■ (1-9) 

1=0 

where Ilij — \i) In the following, we will note oJi^i+i — t^i+i ~ i^i- 

A transmon qubit is a Cooper-pair box operated in the limit where Ej{^) Eq. 



This has important consequences that are discussed extensively in (Koch et al, 2007). 
The transmon is essentially insensitive to charge noise, which increases its coherence 
time considerably. Its eigenfrequencies iOi obey approximate simple relations /iwoi ~ 
y/2Ej{<S>)Ec and ~ — Ec/{4:h). Its wavefunctions are very similar to 

those of a harmonic oscillator. This implies interesting properties for the transmon 
matrix elements: approximately, we have {i\N\j) Sj,i+i\/i -f 1 (0| |1) (assuming 
j > i) and similarly {i\ 9 |j) w 5,^,+l^/T+l (0| 

Since spontaneous decay of the excited states of the transmon occurs mainly by 
linear coupling of the charge or phase degrees of freedom to the various baths (either 
electromagnetic or microscopic), this in turn implies that relaxation of the transmon 

^These relations are approximate. In the numerical simulations presented in the chapter, the exact 
energies and couplings were computed by diagonalizing the Cooper-pair box Hamiltonian. 
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Fig. 1.4 Coupled qubit-nonlinear resonator system. The transmon qubit is capacitively cou- 
pled to one end of the nonlinear resonator, while the other end is coupled to a 50Q, measure- 
ment line through which microwave pulses of amplitude (resp. £3) and frequency Up (resp. 
uus) drive the resonator (resp. the qubit). 

occurs mostly from level i + 1 to level i (and not to lower-lying levels) with a rate 
7i+i_i that verifies 7i+i,i ~ {i + 1)7 where 7 is the spontaneous emission rate of 
level |1) towards |0). Other uncontrolled processes include low- frequency noise (flux 
or charge noise mostly) which causes dephasing of level \i) at a rate 7^^^. All these 
processes are described in the following master equatiorj^ 



M-2 
i=0 



M~2 







(1.10) 



where we considered the qubit to be at zero temperature. In addition, the qubit state 



can be manipulated by driving through the resonator [sec Eq. (1.6)] 



1.1.3 Coupling 

In our setup the transmon qubit is capacitively coupled to the resonator by the gate 



capacitor Cg and located at one end of the resonator as depicted in Fig. 1.4 It is 



straightforward to show (Koch et al., 2007) that such a capacitive transmon-resonator 



interaction is described by a multilevel Jaynes-Cummings model 



M-2 



(1.11) 



with gij = 2e/3 {i\ N \j) SVq, SVq = w^^/iZq/tt the vacuum fluctuations of the voltage 
between the resonator center conductor and ground plane at the qubit location, and 
/3 = Cg/{C + Cg). Note that we have neglected counter-rotating terms because they 
oscillate rapidly and only contribute to a slight energy shift. Due to the properties of 
the transmon matrix elements discussed above gi j ~ + Ig, where we have 

noted g ~ go.i- This yields the following coupling Hamiltonian 



^Note that this description for dephasing is insufficient to describe dephasing of the levels i > f . 
A quantitative model of dephasing for these higher levels should take into account the spectrum of 
the bath causing dephasing. This description is however sufficient for our purpose in this chapter. 



Semiclassical dynamics of the nonlinear resonator 7 



M-2 



A 9^,i+l{a^^^,i+l + an,+i,,). (1.12) 



All the experiments reported in this chapter will be carried out in the limit where the 
detuning A — luqi — lo^ between the transmon and resonator transition frequencies is 
much larger than the coupling constant g, the so-called dispersive limit. Approximate 
expressions for the interaction Hamiltonian can then be derived by applying perturba- 
tion theory to various orders. In this work we will be both interested in a "low-energy" 
limit in which the average number of photons in the resonator n is small compared 
to the critical photon number ricrit = (^)^/(45^), and also in the "high-energy" limit 
where this is no longer true. In this case a more refined treatment is needed (see section 



1.4). In the low-energy limit n ^ ?T.crit, the simplest second-order perturbation theory 



yields an effective Hamiltonian which, restricted in the {|0) , |1)} subspace, writes: 

Hint = Kx<7z + s)a)a, (1-13) 

where cr^ is the usual Pauli matrix (such that the qubit Hamiltonian writes H^ = 
-huimjla^), X = X01-X12/2, xoi = X12 = gl2/{^i2~^r) and s = -X12/2 (|Koch 



et al, 2007[ ). 

In summary, including the dissipative terms, the qubit drive, and the nonlinear 
resonator pump, we get the following master equation 

p = Cp=-'-[Hi,p\+Cqp + CrP (1.14) 
that describes the dispersive interaction between a superconducting qubit and a KNR. 

1.2 Semiclassical dynamics of the nonlinear resonator 

In this section we temporarily leave the qubit aside and study the dynamics of the 



pumped KNR alone described by Eq. (1.7). We identify the bistability region by calcu- 
lating the classical response of the KNR to the pump. We then focus on calculating the 
switching probability of the resonator between its two dynamical states as a function 
of the amplitude and frequency of the microwave pump pulse. Numerical simulations 
are finally compared to measurements. 

1.2.1 Classical response to a pump 

The steady-state complex amplitude a of the intracavity field in the classical limit is 



Yurke and Buks, 2006 ) 



obtained from Eq. (1.7) (Dykman and Krivoglaz, 1979 Drummond and Walls, 1980 



+ -a = -icp, (1.15) 



where 51 — 2(3(1 — Up/ur) is the reduced detuning of the pump and ep its reduced 
amplitude. This equation can be solved numerically, yielding the amplitude of the 
oscillating current / = yJKpKZQU}j.\a\ passing through the Josephson junction. This 



is illustrated in Fig. 1.5 for typical sample parameters. At low drive amplitude e^, 
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Fig. 1.5 Classical response of a KNR to a pump microwave signal, (top) Calculated max- 
imum current I through the junction as a function of the reduced detuning fl for various 
pump powers in critical power (Pc) units, (bottom) Stability regions of the oscillator driven 
at frequency Q and power Pp. Region L (H) corresponds to low (high) oscillation amplitude, 
whereas in the hatched region both solutions are metastable. In the bistable region {Q > Qc) 
the KNR acts as a Josephson Bifurcation Amplifier (JBA), whereas for < Sic the KNR 
acts as a Josephson Parametric Amplifier (JPA) with maximum gain close to the critical 
point (QcPc)- 

the response is a Lorentzian centered around = as for a linear resonator. When 
the drive amplitude increases, the resonance shifts downwards and shows a sharpened 
response for some frequencies. In this regime, small variations in the pump amplitude 
produce a large change in a, indicating that the KNR acts as a Josephson parametric 
amplifier (JPA) ( jCastellanos-Beltran and Lehnert, 2007[ ). For some critical drive power 
Pc, the slope \da/d€p\ becomes infinite at = flc = a/S- For Pp > Pc, two stable 
solutions L and H with different oscillation amplitude a exist. In this bistable regime, 
occurring only for il > Qc = a/S, the transition from L to H occurs abruptly when 
ramping up the pump power at a value P+(r2) called the bifurcation threshold and 



is hysteretic (Landau and Lifshitz, ). This regime is called the Josephson Bifurcation 



Amplifier (JBA) (Siddiqi et al., 2006c|. Fig. 1 1.5 [summarizes these properties. 



1.2.2 Fluctuations: a qualitative discussion 

Fluctuations are essential to understand the physics of a pumped KNR. When the res- 
onator is biased in the JPA regime, fluctuations can be either amplified or de-amplified 
depending on their phase relation with the pump. This results in a phase-dependent 
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noise of the parametric amplifier. When the resonator is biased in the JBA regime, the 
switching from the L to H state is a stochastic process because of fluctuations which 
make the dynamics of the intracavity field differ from one experiment to the other 
even with identical biasing conditions. A quantitative account of fluctuations is then 
needed in order to determine the probability that the resonator switches from L to H, 
as well as the switching rate. 

In the JPA regime, fluctuations can be treated by linearizing the equations of mo- 
tion around the steady-state solution a [see Eq. (1.151]. The input-output relations are 



then solved showing how input fluctuations are amplified or de-amplified at the out- 
put, with a gain that depends on the pumping conditions. This method has been used 



for instance in (Drummond and Walls, 1980 Collett and Walls, 1985 Yurke and Buks, 



2006 1 to derive explicit expressions for the phase-dependent gain and noise properties 



in the JPA regime, and a complete quantum theory was developed in (Dykman et al., 



2011). In the classical version of this problem, such analysis as well as the effect of 



switching events on the gain and on the power spectra was considered in (Dykman 



and Krivoglaz, 1979 Dykman et al., 1994). Calculation of switching rates in the JBA 



regime (which is our main interest here) requires going beyond a linearized analysis 
of the problem. Approximate theoretical expressions for the transition rate (Dykman 



and Krivoglaz, 1979 Dykman and Smelyanskiy, 1988) can be obtained by the so- 



called real-time instanton technique applied to the full non-linear system. Using the 
critical slowdown of the KNR dynamics close to the bifurcation threshold the problem 
can be mapped onto the Fokker-Planck equation describing the escape of a classical 



Brownian particle in a ID cubic potential (Dykman and Krivoglaz, 1980 Drummond, 



1986 Dykman, 2007). At this point one can obtain expressions for the switching rate 



using standard techniques. A very accessible account of this approach is given in ( Vi 



jayaraghavan, 2008). However these approaches only treat the case in which the pump 



of the KNR is kept at a constant amplitude or varied very slowly compared to the 
KNR dynamics. Experimentally, the switching probability is measured with microwave 
pulses that can induce transients of the intra-resonator field when the rate at which 
the pump amplitude is varied is comparable to the ringing time of the KNR. For 



the parameters of our experiments (see sections 1.3 and 1.4), the duration of the pulses 
indeed requires to deal with transients. 

For that reason, we rely here on numerical simulations of the KNR switching from 
L to H when it is subject to a microwave pump pulse. This approach can treat arbi- 
trary dynamics of the pump amplitude and naturally takes transients into account. 
In the following, we compare measurements of the switching curves of a KNR whose 
parameters are known with great accuracy to the simulations. 



1.2.3 Experimental implementation 



The experimental setup is sketched on Fig. 1.6 a). The KNR is capacitively coupled 



to a transmon qubit. In the experiments reported in this section, the qubit is biased 
at a flux such that its resonance frequency is much lower than the resonator's to avoid 
perturbing the KNR dynamics. The KNR has a fundamental frequency uJr/2TT = 6.4535 
GHz (in the absence of qubit), quality factor Q = 685 and Josephson junction critical 
current Iq — 720 nA. These parameters are determined experimentally in separate 
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Fig. 1.6 (a) Experimental setup for S-curve measurement. A microwave pulse at frequency 
ujm is sent onto the KNR. After reflection onto the KNR, the pulse is amplified and de- 
modulated at room-temperature. The oscillogram shows filtered I{t) traces corresponding 
to switching and non switching events for Pm close to the bifurcation. The dashed line is a 
threshold defined to extract the switching probability Ps and acquire S-curves. (b) Temporal 
profile of the microwave pulse: the microwave amplitude is ramped linearly up to power Pm 
in a time tr, kept constant for a measurement time tm, and reduced to power Ph during a 
hold time tu- (c) Typical S-curve Ps{Pm)- 

measurements. Microwave pulses at frequency uJm (corresponding to reduced detuning 
Q) are sent to the KNR with a shape shown in Fig. 1.6 b). The microwave amplitude 



is first ramped linearly in tr from to the maximum amplitude corresponding to a 
power Prn- The power is then kept constant for a time tm (effective measurement 
time) during which the resonator may or may not switch from L to H. The power is 
finally decreased to a hold value Ph at which the resonator state cannot change due to 
hysteresis, for a hold time t^ long enough (typically 1/Us) to allow good discrimination 
at room temperature of L and H (Siddiqi et al., 2006c| ) . After the pulse is reflected by 
the KNR, it is amplified, and the two quadratures / and Q are measured by homodyne 
detection at room temperature. For a power Pm close to the bifurcation threshold, 
the time traces belong to two clearly resolved families of trajectories corresponding 



to oscillator states L and H [see Fig. 1.6 a)]. By defining a measurement threshold 



[dashed line in Fig. 1.6 a)] and repeating the measurement sequence, we can extract 
the probability ps to find the resonator in state iJ as a function of the drive parameters 
Pm and ft. We thus obtain so-called S-curves as displayed in Fig. 1.6 c). 



1.2.4 Comparison to numerical simulations 

Numerical simulations of the S-curves are performed as indicated in the following. The 
master equation ( 1.71 is integrated in the Fock state basis, starting with the resonator 
initialized in the vacuum state (or in a thermal state at temperature T for finite 
temperature calculations), and with a time- dependent drive amplitude ep(t) taken 
in exact correspondence to the experimental pulse profile. At the end of the hold 
time tn, the function Q{a) = (Q;|p|a) {Q distribution (Walls and Milburn, 2008)) 
is computed from the density matrix p. When the pulse amplitude approaches the 
bistability threshold, the Q distribution shows two well-separated regions [see inset of 
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Fig. 1.7 a)] corresponding to resonator states L and H. The probability ps is then 
given by the integrated weight of the high-amphtude part of the Q distribution [see 
inset of Fig. 1.7 a)]. The temperature of the bath can be easily varied as described 



in (1.7). If the bath is at zero temperature, only quantum fluctuations are present. 




Fig. 1.7 Comparison of measured and simulated S-curves. (a) S-curves at T=20 mK for 
different Q. Squares are experimental data, dots are simulated data (the grey line is a guide 
to the eye for the latter). Inset: greyscale representation (black=0, white=0.15) of the function 
Q{a) = {a\ p \a) used to infer the circled simulated point in the center of the main panel, p 
being the density matrix of the intra-resonator field calculated by numerical integration of 
the master equation, (b) Temperature dependence of the S-curves widths for — 8.2. Dots 
linked by a line are obtained by fitting the simulated S-curves with an erf function (see text). 
Squares are obtained by fitting the upper part of the experimental S-curves as shown in the 
inset. 



In Fig. 1.7 a) we first show three S-curves measured for different Q with a slow pulse 
tj. = 300 ns and t^^ — 2 ^s, together with the corresponding simulations rescaled by a 
global attenuation factor that we have adjusted to fit the data. We first note that the 
position of the three S-curves is perfectly reproduced by the simulation, indicating that 
the parameters of the simulation are compatible with our sample. The overall shape 
of the S-curves is also well reproduced by the simulation, despite a structure at low ps 
which is present in the experiment but not in the simulations and that we attribute 
to the residual effect of the transmon qubit on the KNIi[^ Apart from this structure 
at the S-curve foot, the agreement between the simulations and the experiment is 
remarkable at all VL with the global attenuation as the only adjustable parameter. 



^Indeed, in these measurements the transmon was dcUberately biased at low frequency (around or 
below 4 GHz to shift it far away from the KNR frequency), resulting in a finite thermal population 
even at the lowest temperature T = 20 mK of our cryostat. Excited states of the transmon correspond 
to a S-curve shifted at a lower Pm , and thus the observed S-curve should be seen as a sum of several S- 
curves weighted by the various excited state populations causing the structured S-curve foot observed 
in Fig. |1.7| a). We are unfortunately unable to quantitatively account for this structure since in this 
experiment we don't accurately know either the transmon's frequency or its effective temperature. 
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Since these simulations were realized assuming zero temperature, this constitutes a 
strong indication that the width of the experimental S-curve is essentially dominated 
by quantum fluctuations. Note that we expect that the transitions between the two 
states are not due to tunneling but take place over the effective potential barrier, in 



the process termed 'quantum activation' in ( Marthaler and Dykman, 2006 ) 



We also measure the temperature dependence of the S-curve width (determined 
taking into account only the large-ps part of the curve) for similar pulse parameters 
and compare it to the simulations [see Fig. 1.7 b)]. The agreement is quantitative up to 
T K, 200 mK. Deviations are observed between 200 and 300 mK which might be caused 
by a decrease of the KNR quality factor due to thermally excited quasiparticles, or by 
thermal excitation of the qubit. We observe that the crossover from the thermal to the 
quantum regime in both the experiment and the simulations occurs around 100 mK. 
This is significantly lower than the theoretical prediction of huJr/2kB which is around 



150 mK for our parameters (Marthaler and Dykman, 2006 Vijay et al., 2009). More 



detailed investigations are needed to address this question. 



1.3 High-fidelity qubit readout by a nonlinear resonator 

We now describe how the bistability of the nonlinear resonator combined with its 



dispersive interaction with the qubit [see Eq. (1.13l] can be exploited to yield a high- 
fidelity qubit state readout ( Mallet et al., 2009| ). 



1.3.1 Principle 

The underlying principle of this readout method can be understood with the dispersive 
qubit-resonator interaction Hamiltonian (1.131, valid provided |A| and n ricrit 



as explained in the first section. In this regime, the resonator frequency takes the qubit- 
state dependent value ujr^i = ujr + s ± x- Detecting this frequency change sufficiently 
fast before the qubit relaxes by spontaneous emission is the key challenge of all readout 
methods in circuit QED. This is usually done by sending a microwave pulse close to 
the resonance frequency and measuring the phase of the signal transmitted through or 



reflected by the resonator (Blais et al., 2004), which bears information on the resonator 



frequency and therefore on the qubit state. One outstanding feature of this dispersive 
scheme is that it perturbs the qubit state minimally as evidenced by the fact that 
the interaction Hamiltonian commutes with az, implying that the readout should be 



quantum non-demolition (QND) (Braginsky and Khalili, 1996). 

However this readout method suffers from a limited signal-to-noise ratio (SNR) 
when applied to a simple linear resonator. Indeed the mean photon number n in the 
resonator should be kept much lower than ricrit = A^ / (4g^ ) to avoid spurious qubit 



excitation ([Boissonneault, Gambetta and Blais, 2008[ Boissonneault, Gambetta and 



Blais, 2009), which limits the measurement power. Such a small signal must be am- 



plified for detection but state-of-the-art cryogenic amplifiers have a noise temperature 
Tn « 3 K so that the added noise is much larger than the signal. Typical figures are 
g/A ~ 1/10, so that ricrit ~ 25 whereas the amphfier adds a noise corresponding to 25 
photons. Averaging during a time improves the SNR, but the longer tm, the higher 
the probability that the qubit (characterized by a decay time Ti of the order of the 
microsecond) relaxes yielding an incorrect readout. There thus exists an optimal t^^ 
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yielding an optimum discrimination between both qubit states which can be shown to 
be of order Ti (Gambetta et al., 2007). These constraints have up until now made it 
impossible to reach the readout efficiencies required for quantum information applica- 
tions, even though quantum state tomography is still possible by performing ensemble 
averages (Filipp et al., 20091. In the following, we show that turning the resonator 
into a KNR makes it possible to improve the visibility of the measurement without 
compromising either its QND character or the qubit coherence properties. 




Fig. 1.8 Principle of JBA readout, (a) Sketch of the stability diagrams of the NL resonator 
dispersively coupled to a transmon prepared in state |0) (black) or |1) (grey).(b) Sketch of 
the probabihty ps to measure the resonator in state H for qubit initiaUy in |0) (black) or |1) 
(grey). 



Taking advantage of the bistability, a KNR dispersively coupled to a qubit can 



be turned into a sample-and-hold qubit state detector with a large SNR (Siddiqi 



et al., 2006a Boaknin et al., 2007). Due to the dispersive coupling (1.13) the re- 
duced pump detuning fli and the bifurcation threshold P+(rJi) now depend on the 
qubit state. Sending microwave pulses at a frequency and a power Pm such that 
P+(r2i) < Pm < P+(rio) then maps the resonator states L or H onto qubit states 
|0) or |1) as shown in Fig. 1.8 a) (Mallet et al., 2009). The pulses have a temporal 



shape identical to those described in the previous section [see Fig. 1.6 b)], allowing 
perfect single-shot discrimination of the two resonator states as already discussed. 
The key advantage of this readout method over readout with a linear resonator is 
that the outcome of the measurement is totally insensitive to relaxation of the qubit if 
the qubit relaxes during the hold part of the readout pulse. However qubit relaxation 
occurring before the resonator switches causes measurement errors as we will see be- 
low. This approach results in potentially large readout fidelity provided the separation 
between the S-curves corresponding to both qubit states is large compared to their 
width [see Fig 1.8 b)]. Such an rf-switching scheme has been implemented successfully 



to readout a quantronium qubit (Boulant et al., 2007 Siddiqi et al., 200661, a flux 



qubit (Lupa§cu et al., 20061 and a transmon qubit (Mallet et al., 20091, as described 
in the next section. We finally note that two recent experiments use nonlinearity to 
achieve high-fidelity qubit readout in a different way: Vijay et al. (Vijay et al., 2011) 
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have used a JPA as a first-stage amplifier for linear dispersive readout and in this way 
have observed quantum jumps of a transmon qubit, while Reed et al. (Reed et al, 



2010 1 use the nonlinearity of the cavity inherited by coupling to the qubit to obtain 



qubit-state dependent signal in the ultra-strong driving regime ( Boissonneault et al., 



2010 Bishop et al., 2010) 



1.3.2 Implementation of the JBA readout of a transmon 
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Fig. 1.9 Readout contrast: experiment and simulation. S-curves for qubit prepared in |0) 
(black) or |1) (dark grey) with ujm/'in = 6.437 GHz = aj„ii/27r. Full lines are data and dashed 
lines are simulations. The black arrow shows the maximum experimental contrast. The lighter 
grey full line is the S-curve obtained with qubit prepared in |0) with the same parameters as 
for the solid lines except that uurn = i^mi — Aoji (see text). The grey arrow shows the expected 
contrast in absence of relaxation. 




The experimental setup is the same as shown in Figs. 1.4 and 1.6 In contrast to 
the measurements reported in section [1.21 the transmon qubit is now biased at a flux 
such that its frequency is closer to the KNR so that the dispersive coupling constant 
X becomes sizeable, a necessary requirement for readout with high contrast. In all the 
measurements reported here the qubit is operated at negative detunings A — cjoi — Wr, 
larger in absolute value than 5/27r = 44 MHz. Readout pulses [Fig. 1.6 b)] of frequency 
LJm and power are sent to the KNR, allowing to measure S-curves as explained 
in [T231 

Figure[r9|shows S-curves obtained at -A/2n = 380 MHz with a;„/27r = 6.437 GHz ; 
ujmi/2TT for both qubit states. The |1) state was prepared with a tt pulse applied just 
before the readout. The maximum separation, or best contrast, at that working point 
is 88 %. Also shown is a S-curve with the qubit prepared in |0) and measured in the 
same conditions but with ujm = i^mi — Awi. The detuning Aloi/2t: = 4.1 MHz is cho- 
sen such that the corresponding S-curve matches best the S-curve with qubit prepared 
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in |0) and w,„ = (jJmi- Thus, Awi yields a direct measurement of the cavity puh and 
is in good agreement with the theoretical value of 2x/27r = 4.35 MHz calculated for 



these parameters (see 1.1.31 



However comparing the grey and the light grey S-curves in Fig. 1.9 we note a 
discrepancy when ps approaches unity. This is due to relaxation of the qubit state 
|1) during the time the resonator energizes and choses its state. Quantitative analysis 
(not presented here) shows that this time is of the order of Ti/10 ~ 40 ns with our 
circuit parameters. Although the consequences are much less dramatic than in linear 
circuit-QED, the JBA readout fidelity remains limited by finite Ti. To emphasize this 
point we estimate that the contrast without relaxation would be 97% (gray arrow on 



Fig. 1.91 



We also present in Fig. |1.9| the results of a simulation of the complete master 



equation (1.14), taking into account the measured value of the qubit relaxation time 
Ti and of the qubit-resonator coupling constant. Even though the shift of the S-curves 
between states |0) and |1) is slightly less in the simulation than in the experiment, the 
overall agreement with the experimental data is remarkable with the global attenuation 
of the measurement line as the only adjustable parameter. In particular the maximum 
contrast predicted by the simulation is precisely the one found in the experiment. This 
demonstrates that despite the complexity of our circuit (a multilevel system coupled 
to a KNR) it is possible to obtain a quantitative and detailed description. 

In the next subsection we show how we can partially get rid of the relaxation to 
improve the fidelity of the readout, and we finally present a general figure of merit of 
the JBA transmon readout. 

1.3.3 Performance and limitations of the transmon readout by a JBA 

To reduce the effect of relaxation and improve the readout contrast we use a technique 
analogous to electron shelving in atomic physics and already used with other Josephson 



qubits (Martinis et al., 2002) : we transfer state |1) into the next excited state |2) 
with a resonant 7r-pulse just before the readout pulse. This has two effects on the 
readout performance. First, it can be shown that the resonator frequency is slightly 
more shifted by state |2) of the transmon than by state |1), yielding a larger S-curve 
separation between the two qubit states and thus an enhanced readout fidelity. Second, 
relaxation from state |2) occurs only towards state |1) (which hardly changes the 
outcome of the readout if Pm is properly chosen), and not directly towards |0) in 



the transmon, as explained in section [1.1.2 Still working at — A/27r — 380 MHz 



we obtain the solid grey S-curve on Fig. 1.10 a) yielding a 92% contrast. The right 



panel shows Rabi oscillations between |0) and |1) obtained with such a composite 
readout pulse. This contrast, larger than 90%, is in agreement with the width of the 
S-curves estimated from numerical simulations, with their theoretical displacement, 
and with the measured qubit relaxation time. Of the remaining 8% readout error, we 
estimate that about 6% are due to relaxation before switching and 2% to residual 
out-of-equilibrium population of |1) and to control pulse imperfections. 

We now discuss the dependence of the readout contrast and qubit coherence on 
the detuning A = wqi — i^r- Besides acting as a qubit state detector, the resonator 
also serves as a filter protecting the qubit against spontaneous emission into the 50 
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Fig. 1.10 Performances of JBA readout, a) Best single siiot readout contrast. Left: S-curves 
witli qubit prepared in jO) (black), jl) (dotted grey) and j2) (solid grey). Right: Rabi oscilla- 
tions measured with the composite readout (dots). The solid line is a fit by an exponentially 
damped sine curve (peak-to-peak amplitude 94% and decay time 0.5 ^s). b) Trade-off be- 
tween qubit coherence and readout fidelity. Top: experimental relaxation time T\ (dots) and 
dephasing time (open circles) of the qubit as a function of ojoi . The solid line is the value 
of Ti obtained by adding to the expected spontaneous emission through the resonator a re- 
laxation channel of unknown origin with T\ = 0.7 /is. The dashed line is the pure dephasing 
time corresponding to a l/oj flux noise with an amplitude set to 20 /i<l?o/\/Hz at 1 Hz. 
Bottom, left axis: Readout contrast with (dots) and without (open circles) transfer from state 
|1) to 1 2). Bottom, right axis: measured (open squares) and predicted (dashed line) cavity 
pull 2x- The grey area highlights the detuning region where the readout contrast is higher 
than 85 %. 



impedance of the external circuit (Esteve et al., 1986 Houck et al., 20071. The smaller 
|A|, the stronger the dispersive coupling x between the qubit and the resonator, im- 
plying a larger separation between S-curves curves but also a faster relaxation due to 
the Purcell effect. We thus expect the contrast to be limited by relaxation at small 
|A|, by the poor separation between the S-curves at large |A|, and to exhibit a maxi- 
mum in between. Figure 1.10 b) presents a summary of our measurements of contrast 
and coherence times. At small |A|, Ti is in quantitative agreement with calculations 
of the spontaneous emission through the resonator. However it presents a saturation, 
similarly as observed in several circuit-QED experiments (Houck et al, 2008). The 
cavity pull determined from the S-curve shifts [cf. Fig. 



1.9 b)] is in quantitative agree- 



ment with the value of 2% calculated from the sample parameters. The contrast varies 
with A as anticipated and shows a maximum of 92% at — A/27r = 380 MHz, where 
Ti = 0.5 /iS. Larger Ti can be obtained at the expense of a lower contrast, and recipro- 
cally. Another important figure of merit is the pure dephasing time which controls 
the lifetime of a superposition of qubit states (Ithier et al., 2005 1. is extracted from 
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Ramsey experiments and shows a smooth dependence on the qubit frequency, in qual- 
itative agreement with the dephasing time deduced from a l/w flux noise of sp ectral 
density set to 20/x<1>o/-\/Hz at 1 Hz, a value similar to those reported elsewhere (Well- 



stood et al., 1987). To summarize our circuit performances, we obtained a 400 MHz 



frequency range [grey area on Fig. 1.10 b)] where the readout contrast is larger than 
85%, Ti larger than 300 ns, and larger than 700 ns. Further optimization of the 
JBA parameters Iq and Q might increase this high-fidelity readout frequency window. 



We finally discuss the QND character (or projectiveness) of the JBA readout (Bra- 



ginsky and Khalili, 1996). Since the dispersive interaction commutes with the qubit 



Hamiltonian, the qubit initially prepared in an arbitrary superposition should be 
projected onto the state that is measured. We define the QND fidelity as Fqnd = 
(P(0|0) + P(l|l))/2, where P(i|i) {i = 0, 1) is the probability th at the qubit state \i) 
is unchanged after a single measurement (Lupascu et al., 20071. Josephson Bifurca- 
tion Amplifier implemented to readout a flux qubit demonstrated a QND fidelity of 



(Lupascu et al., 2007). In the case of a transmon, the data presented in (Mallet 
et al., 20091) for two successive readouts yield a QND fidelity of ~ 40%. However this 



poor QND character was shown to be only due to relaxation occuring during the first 
readout pulse, and not to a lack of projectivity of the measurement. This is due to the 
rather long duration of the measurement sequence [Fig. 1.8 b)]: ts (th) is typically a 



few tens (hundreds) of nanoseconds so that the overall readout time is a significant 
fraction of Ti. As a comparison, the readout pulse duration in (Lupascu et al., 2007) 
was a few tens of nanoseconds, which was possible because of the low quality factor 
{Q w 20) of the resonator. Using such a low Q cavity for the JBA would certainly help 
performing the readout faster but at the price of degrading the fidelity of the readout 
(broadening of the S-curves) and the coherence times (inefficient Purcell protection). 



1.4 Backaction of a driven Kerr resonator on a qubit 

In the readout method described in the previous section, the KNR is biased very close 
to its bifurcation threshold. The small frequency shift due to the qubit is then sufficient 
to push the resonator above the threshold, resulting in a very large difference at the 
output for the two qubit states. In other words, the nonlinear resonator energized by 
the pump signal behaves as an active measurement device with some power gain: an 
on-chip amplifier. This is in strong contrast with the more traditional measurement 
process in linear circuit QED where the resonator has a purely passive behavior. It 
then seems very interesting to investigate the quantum backaction of this specific 
measurement process on the qubit and to compare it to the well-known linear circuit 
QED case. 

Indeed, any measurement in quantum mechanics necessarily induces decoherence 
in the variable conjugate to the one being measured. Quantum physics imposes the 
inequality 

> r„ieas/2 (1-16) 



between the system's measurement-induced dephasing rate and the measurement rate 
Lmeas, Stating that the most efficient detector can only measure as fast as it de- 
phases (Clerk et al., 2003). The dispersive readout in cQED with a passive linear 
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resonator has been shown theoreticaUy to reach this quantum hmit (Gambetta et al., 
[2008). In the fohowing, we will investigate whether the qubit readout by an active 
nonlinear resonator reaches that limit. We will first give a theoretical treatment of 
this question and then present an experiment performed on the same sample as used 
in the previous section. In this experiment, we have measured the qubit dephasing 
rate as a function of the drive amplitude and frequency of the nonlinear resonator, 
allowing us to test these predictions. 

1.4.1 Backaction of a nonlinear resonator on a qubit: theory 

Backaction in linear cQED: theory. Before taking into account the Kerr nonlinearity 
of the resonator, we first review the linear case and consider only two transmon levels 
for simplicity. Let us consider the qubit prepared in the superposition cq |0) + ci |1) 
; one then turns on a measurement drive at a frequency Wp chosen to be close to 
the resonator frequency w^. The dispersive interaction Eq. (1.13) then leads to an 
entangled qubit-field state of the form 



(1.17) 



with Ci_i{t) ~ |cjp as required for a QND interaction. In steady-state, the qubit-state 
dependent coherent state amplitudes are given by 



"0,1 



{uJr — LOpztx) — iK/2 



(1.18) 



and are illustrated in Fig. 1.11 For large enough distances D = \ao — ai\, these two 



field states can be distinguished in a homodyne measurement and thus act as pointer 
states for the qubit. The parameter D is called distinguishability and plays a key role 
in the following. 

From quantum trajectory theory, the optimal rate at which this measurement ex- 
tracts information about the qubit state is found to be T^eas = i^D^ and is propor- 
tional to the intracavity photon number n (Gambetta et al., 2008). Moreover, the fun- 



damental backaction of the measurement on the qubit is found in the time-dependence 
of the off-diagonal elements of p(t), which, neglecting transients, decay at a rate 



r^m = 2xIm[aoai] = -D^ 



(1.19) 



saturating the inequality Eq. (1.16). This expression was experimentally confirmed by 



spectroscopic measurements of the qubit linewidth as a function of n (Schuster et al., 



2005[ |Palacios-Laloy et al, 2010| ). 



The off-diagonal elements of p(t) also experience a frequency shift 6oJa = 2xRe[ao{t)al (t)], 
corresponding to the ac-Stark shift. In the limit 2% < k, this shift linearly follows the 
intra-cavity photon number Suia ~ 2xn (Schuster et al., 2005). The increase of the 



qubit linewidth with intracavity photon number can then simply be understood 
as photon shot-noise in the ac-Stark shift. Indeed, each number state \n) entering the 
pointer states ja^) is shifting the qubit frequency by a quantity 2xn resulting in a 
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Fig. 1.11 Graphical representation of tlie pointer states \ai} in phase space and of the 
distinguishability D = |ao — aij. 



broadened spectroscopic line. When 2x > k, the photon number states can be re- 



solved spectroscopically leading to number splitting ( 


Dykman and Krivoglaz, 1987 


Gambetta et al, 2006 


Schuster et al, 2007 


Clerk and Utami, 2007 


I. 



We note that the classical fields ao(i) in this linear case can be obtained by lin- 
earizing the response around an average field a calculated without the qubit signal 
provided 2% <^ k. Such a linear response theory can still hold for a nonlinear resonator 



at low pump power and this regime was studied before (Laflamme and Clerk, 2011). 
Here, we are interested in regimes in which the bias (pump) drive can approach the 
bifurcation threshold. Close to these points, the susceptibility Ida/dSlj — and hence 
the gain \da/dep\ — of the KNR diverges. The range of signal amplitude (the qubit 
cavity-pull) for which linear response theory holds therefore vanishes. 

Backaction in nonlinear cQED: theory. We now turn to the backaction of the non- 
linear resonator on the transmon qubit, taking into account its multilevel structure. 
To go beyond linear response, we use the polaron transformation approach developed 



in (Gambetta et al, 20081. This approach was first used for a two-level system dis- 
persively coupled to a linear resonator. To be useful with a nonlinear resonator and a 
many- level qubit, it must however be adapted in several ways. 

First, in the dispersive readout presented in section [T. 4. 1[ the photons inside the 
resonator are assumed to have a frequency w^. This holds for a linear resonator because 
one always measures with a drive that is (quasi-)resonant with ujr. On the contrary, 
we saw in section |1.3| that the JBA readout drive is usually detuned from ojr by a 
few tens of megahertz, or equivalently a few percent of the qubit-resonator detuning. 
This detuning yields a variation on the ac-Stark shift per photon that needs to be 
accounted for to obtain quantitative agreement between the theory and the experi- 
mental data. Meanwhile, the qubit is still Lamb-shifted by vacuum fluctuations at the 
resonator frequency. This frequency however now depends on the amplitude of the 
field inside the resonator. Moreover, we are interested in the qubit spectroscopy while 
the resonator being pumped. We must therefore consider two drives on the system: 
one far-detuned from (pump) and one quasi-resonant to (spectroscopy) the qubit fre- 
quency. Finally, and as stated above, we must go beyond linear response theory. In the 
following subsections, we develop a theory ( Boissonneault et al., 20lT] ) that accounts 
for the aforementioned effects. 
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Classical field. The first and most important effect to consider is tlie classical field 



created by the pump drive. A displacement transformation in the bare frame (Blais 



et al., 2007) and a polaron transformation in the dispersive frame (Gambetta et al., 
2008 1 are two approaches that have previously been used to treat this field. The former 
captures the pump-resonator detuning in the ac-Stark shift, but, since there is only one 
field a, it is within linear response theory. On the contrary, the latter can go beyond 
linear response but, since the qubit-resonator dispersive approximation is done first, 
the ac-Stark shift is a function of the qubit-resonator detuning rather than of the 
qubit-pump detuning. 

Here, we first do a polaron transformation on the full Jaynes-Cummings Hamilto- 
nian. The polaron transformation is defined as 

M-l 

P=Y1 n,,,D[a,(0], (1.20) 

i=0 



where D[Q!i(i)] is a time-dependent field displacement operator (Gambetta et al., 



2008). Defining H' = P^HP, we obtain the Hamiltonian in the polaron basis 

H' = H'j^L + H'a + H'g + H'j + Hp = H'o+ H[+ H'^, (1.21) 

where Hp = iP^P = {—iHaa^ + ill* a) — ImpcII* ] accounts for the time-dependence 
of the polaron transformation with: 

i?o = + I] 3^ [n:n,,,+i + n,+i,,;n„] + s(n„), (1.22a) 

i=0 

M-2 

H[^Y. ft(n,,.+ia^ + h.c.) + [/(n„)at + h.c] , (1.22b) 

i/^ =cj;(n„)a^a, (1.22c) 

where w'^ljla) = + 2i4r|nQp -f S^T'IIIqI^, silia) and /(IIq) are functions that are 
diagonal in the qubit basis and act as identity in the qubit subspace if — aj, and 
IIq, = ^f^Q^ ^i,iCti{t). The subscript on these three transformed Hamiltonians refer 
to the number of resonator ladder operators involved. In obtaining these transformed 
Hamiltonians, we have assumed that the "classical part" of the fields is larger than 
the "quantum part", i.e. that la^p 3> (a^a^_^_^, where this last mean value is taken 
in the steady-state and in the polaron frame. Doing so, we have dropped all terms 
with more than one resonator ladder operator (a^ ja^ja^ aa,...), except diagonal one s 



(a' a). These terms would lead to squeezing-related effects (Laflamme and Clerk, 2011 ), 



which we consider to be small here. As shown in [Gambetta et al. 2008 the polaron 
transformation applied on H, yields 

H', ,,+1 = H,,,+iDt[a,]D[a,+i] = H,,,+iD[a,+i - a,]e-*i-["*+i"-l (1.23) 

Here, we truncate the series at the zeroth order and rather take = H^^^+i. While 

we allow for the pointer states to be different, this simplification imposes the limit 
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a,|2«l/(ata)^^ (1.24) 



on their separation. We note that this approximation can hold for a large separation of 
the pointer states if (a^a) is small in the polaron frame. It is exact in the dispersive 
regime for a linear resonator and also is a good approximation in the limit of small 
squeezing. 

Following |Gambetta et al. 2008| applying the polaron transformation on the dissi- 
pative part of the master equation yields 

V[ay = V[a]p' + 2?[n„] - ^ [i(n* a - n„a^), p'] + a [p\ H^] + [n„, p'] a). (1.25) 

The second term of this equation contains measurement induced dephasing, while the 
third acts as a drive on the resonator and will be cancelled by a proper choice of the 
cki's. The two last terms act as lowering operators in the polaron frame and can be 
neglected if (a) aj^^ is small. In the following steps, we want to diagonalize H'q — which 
will give us the ac-Stark and Lamb shifts — and choose the displacement operator P 
so as to eliminate eliminate H'^ , which will ensure that the photon population is small 
in the polaron frame. 

Stark shift. The effective drive on the qubit giTlaTli^iJ^-i + h.c. yields an ac-Stark 
shift of the qubit frequencies. To compute this shift, we first assume that the field in 
the cavity can be written as IIq, = IIq e~'"p* + IlQ.^e"*'^'*. Since u)p is far-detuned 
from the qubit frequency, this drive can be treated dispersively. However, since uig 
can be resonant with the qubit, this drive must be treated exactly even if is small 
in practice. Here, we obtain the effect of the pump drive on the qubit by applying a 
dispersive transformation of the form 



D = cxp 



'M-2 



J2 6Wn.,^+i"rWn,+i,. 



i=0 



(1.26) 



on the Hamiltonian and on the dissipative Linblad terms of the master equation. 
Considering only the pump drive (since Hq,^ is assumed to be small), choosing 

= , (1.27) 

and expanding H'^ = D^HqD to fourth order using Hausdorff 's relations, we can show 
that the off-diagonal terms related to the pump drive in Hq vanish. We obtain 

H" = H[; + H'{ + H'^ + Fd, (1.28) 



where 
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M-2 



Hr. + Hri — n^" 



(") 



XI *"=*n,+i + h.c.), 



(1.29a) 



M~2 



i=0 

(1.29b) 

\a)^u;, + Sf\a,f + lKf\a,\\ (1.29c) 

5f = -(xf-xLi), (1.29d) 
= -4^f(|ArP + |Ar_in (1.29e) 



and where A? 



f/iAf, and H'r^ 



contrary to (Gambetta et al., 2008) and in agreement with (Blais et al., 2007), the drive 



H!^. We note that 



frequency ujp appears in the dispersive shifts rather than the resonator frequency. This 
is the correct dependence, and resuhs because we made the polaron transformation 
before the dispersive approximation. 

When obtaining Hq and H2, we have assumed that ap,i ~ apj = Up. Still, contrary 
to what is obtained by a single displacement transformation rather than the polaron 
transformation, the field depends explicitly on the qubit state. Indeed, ap.i and Us^i 
are given by 



= fp.,{a) = [{ujr - + 5, - ifj + {K+ l,Ki) |ap,,p + i^Vp.^l^] ap,» + ep, 

(1.30a) 

+ e,. (1.30b) 



= /,^,(a) = [{u 



tf) +if|ap,,f + if'|ap,,|4] 



These choices ensure that the second term of H[ is zero and therefore eliminates the 
photon population in the polaron frame. We neglected frequency mixing terms and 
|app to obtain these conditions. We note that, from the point 



assumed that \a\ 

of view of a signal parametrically amplified by the resonator, Eq. (1.30a) contains the 
qubit signal Si + ;|r^i|Q^p.iP directly as a parameter in the nonlinear equation for ap^i. 
This places this theory beyond linear response. 



Lamb shift. When the conditions Eq. (1.30) are met, the only remaining off-diagonal 



terms in H" beside the spectroscopy drive is the Jaynes-Cummings coupling in 
The only difference between H" and the initial Hamiltonian ignoring the pump drive 
is that the qubit and resonator frequencies now depend on the amplitude of the cavity 
field ttp, through the ac-Stark shift or the Kerr nonlinearity. Since this renormalized 
resonator is not driven, the photon population in this frame is zero and the Jaynes- 
Cummings coupling yields only a Lamb shift. Indeed, diagonalizing using the usual 



dispersive transformation (Carbonaro et al., 1979) and projecting on the two-level 
{|0) , |1)} subspace yields 
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H'" = 



2 

r'" „"'l 



(Tz + go {ctsB *'^^"*(T+ + h.c.) + ujl.{a)a'a, 



^ oz ^ yo {'J^sf^ " <-'+ -I- ii-^-j 

(J" = -^ [H'\ p'"] + i^V[aAp"' + l^-Dla+y + '^V[az]p 



KD[a]p" 



(1.31a) 
(1.31b) 



where 
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Ar(«) 
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c.,:Vi(a)-<(a)-<(a)' 



7; =7 + 7T + |Ao("p)l 

7t = (27^ + «^i?^)|Agnap|^ 



7,^ =7¥> + 



2 2.92 



(1.32a) 
(1.32b) 
(1.32c) 

(1.32d) 
(1.32e) 

(1.32f) 



D ~ |op,i — ctp,o\ is the distinguishabihty between the pointer states of the nonhnear 
resonator. The resuhing master equation represents a qubit for which the hnear cx Sf 
and quadratic cx Kf ac-Stark shifts depend on the pump frequency Wp and the Lamb 
shift LI depends on the detuning between the Kerr-shifted resonator frequency (a) 
and the Stark-shifted qubit frequency. This qubit is subjected to modified relaxation 
rates. 

The rate 7^ contains dressed-dephasing as well as a very similar effect which we 
will refer to as dressed measurement-induced-dephasing. Dressed-dephasing represents 
population mixing caused by the qubit 's pure dephasing bath when the qubit is dressed 
by photons (Boissonneault et al., 2008 ). Dressed measurement- induced dephasing is 
the same effect when the dephasing occurs because of measurement-induced dephasing. 



As shown in (Boissonneault et al., 2009), the precise value of these rates will depend 



on the baths' spectra. The third term of 7^ is Purcell relaxation (Purcell, 1946 Houck 
et al., 2008). In this nonlinear system however, the value of the Purcell rate depends 



on the Kerr-shifted resonator frequency and the Stark-shifted qubit frequency. 



The second term of 7"' is measurement- induced dephasing (Schuster et ai, 2005 



Gambetta et al, 2006 



Gambetta et al., 2008|). We recover the same expression as 

are the solution of the nonlinear 



Gambetta et al. 20081 Here however, the fields a 



equation (1.30a I, in which the cavity-pull Sf depends on qubit-pump detuning instead 



of the qubit-resonator detuning. Finally, the last term of 7^" is dressed-decay (^Boisson- 
neault et al., 2008), and is negligible given the parameters of the experiment described 



in the next section. 



Discussion. Neglecting the last term in Eq. ( 1.32f I which vanishes in the limit of neg- 
ligible spontaneous emission rate 7 yields 7^' = 7,^-1- 



the exact same equation as 

for a linear resonator [see Eq. (1.19)]. This remarkable non-trivial result indicates that 
the direct link between measurement-induced decoherence rate and distinguishabihty 
pertains even when the resonator is nonlinear, provided the approximations made in 
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the previous paragraph are satisfied. Nevertheless the raeasurement-induced dephas- 
ing of a KNR is quafitatively different from the hnear resonator case as we will see in 
the next section. The reason is that the pointer states amplitude, and therefore their 
distance in phase space, have a nonlinear dependence on the pump amplitude and 



frequency as expressed by Eq. (1.30a 



In order to verify if the quantum limit Eq. (1.16 1 is reached here, must be 



compared to the rate Tmeas at which information about the atom state is acquired. 
If the two pointer states are coherent states, then Tmeas — nD"^ as discussed in the 



beginning of this section. Eq. ( 1.32f ) then implies that the backaction of the nonlinear 



resonator is at the quantum limit. According to the quantum theory of the pumped 
KNR however ( [Laflamme and Clerk, 2011 Drummond and Walls, 1980), this is true 
only in the limit of small parametric gain G ~ 1, corresponding to Pp being significantly 
different from P+{il). Indeed, when G ^ 1, the intracavity field shows enhanced 
phase-dependent fluctuations and the field refiected on the cavity shows some degree 
of squeezing. In that limit the measurement rate needs to be re-evaluated and our 
results do not allow us to conclude on the quantum-limited character of the pumped 
nonlinear resonator backaction. 



In related work (Laflamme and Clerk, 2011) Laflamme and Clerk make a different 



set of approximations that allows them to investigate the large gain limit. They assume 
that the qubit-resonator interaction can be described within linear response theory, 
which implies that D is developed to first order in x and that entanglement between 
the qubit and the resonator is neglected. In that framework, they are able to take 
into account the phase-dependent fluctuations of the resonator field and to treat the 
case of arbitrary parametric gain. With these approximations, they show that the 



measurement-induced dephasing rate is still given by Eq. (1.32f). Moreover, they find 
that the measurement rate is Tmeas — kD'^/G because of squeezing in the output field. 
The nonlinear resonator backaction on the qubit thus misses the quantum limit by a 
large factor of order G. 

We note however that typical experimental parameters in circuit QED are usually 
beyond the linear response theory due to the strong qubit-resonator coupling. Fur- 
thermore, the degree of validity of the linear response theory shrinks as the large-gain 



limit is approached, so that the conclusions in (Laflamme and Clerk, 2011) are strictly 



speaking valid only in the limit where x tends towards zero. To conclude, more the- 
oretical work combining a polaron-type approach as in the previous section and an 



account of squeezing such as performed in (Laflamme and Clerk, 2011) seems needed 



in order to assess rigorously the question of the quantum limit of the pumped KNR 
backaction on a strongly coupled qubit. 



1.4.2 Backaction of a nonlinear resonator on a qubit: measurements 



We now present experimental data (Ong et ai, 2011) aiming to test these predictions 



regarding the pumped nonlinear resonator backaction on a qubit that can be summa- 



rized as: 1) a qubit frequency shift given by Eqs. ( 1.32a) and ( 1.30a) 2) dephasing with 



a rate given by Eq. ( 1.32f ). For that purpose, we perform spectroscopy of the artificial 



atom in the presence of a pump field of variable power and frequency driving the KNR. 
The setup, illustrated in Fig. 1.12 a), is similar to the ones used in the previous section. 
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Fig. 1.12 (a) Experimental setup: a transmon qubit is strongly coupled to a coplanar res- 
onator made nonlinear with a Josephson junction. The sample is cooled to 20 mK and is 
driven through an attenuator A by three microwave sources. Source p is used to establish a 
pump field in the resonator, source s for qubit spectroscopy, and source m as a JBA readout: 
its signal at Um is reflected from the resonator and routed through circulators to a cryogenic 
amplifier, a demodulator, and a digitizer, which yields the JBA switching probability ps and 
thus the probability of the qubit excited state (b) Stability diagram of the resonator in 
the Q-Pp plane. The solid lines indicate the highest PA gain below fic and the power at which 
the resonator bifurcates from the low- (L) to the high- (H) amplitude state, respectively. The 
hatched area is the bistability region. Vertical dotted lines correspond the two datasets of 
Figs. [TTsl and [ITil 



In all that follows, the qubit-resonator detuning is fixed at A/2tt = 732 MHz ^ g/^ir, 
so that their interaction is well described by the dispersive Hamiltonian Eq. (1.131 
with x/2Tr — —0.8 MHz and s/27r = 1.7 MHz. Qubit-state readout is performed as 
described in section 1.3 at frequency a;,„/27r = 6.439 GHz {ilo/ftc = 2) by measuring 
the JBA switching probability ps . 

The pulse sequence is shown in the inset of Fig. |1.13[ we apply the spectroscopy 
pulse of varying frequency Ws and fixed power P^, while the resonator is driven by a 
microwave pump pulse of varying frequency ujp and power Pp, followed at the end with 
a qubit readout pulse. Note that the pump pulse starts long before the spectroscopy 
pulse so that the intracavity field has reached its stationary state. Moreover, the read- 
out pulse is applied 200 ns after switching off both pump and spectroscopy pulses. This 
time is long enough to let the intraresonator field relax before readout, but shorter 
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Pp (dBm) 




Fig. 1.13 2D plots of the switching probabihty ps as a function of ujs and Pp for Slo/f^c = 3.1 
and Q,o/Q.c — 0.7 (left and right panels). A few qubit lines are shown in overlay for negligible 
field amplitude in the resonator (top left), near the switching point at Pp = 1.0 dBm for 
Oo/Oc = 3.1, and at Pp = -12.4 dBm and -8.6 dBm for Qq/^c = 0.7. The horizontal 
dashed line indicates the qubit frequency at zero cavity field. Lorentzian fits of the qubit lines 
(see example at —8.6 dBm) yield the ac-stark shift Aa;oi and the FWHM linewidth w. Inset: 
microwave pulse sequence used. 



than the qubit relaxation time Ti = 700 ns. Fig. |1. 13] shows the resulting qubit spec- 
trum as a function of Pp at two pump frequencies above and below flc, respectively in 
the parametric amplifier and the bistability regime. We observe the ac-Stark shift of 
the qubit resonance frequency towards lower values, as well as the broadening of the 
qubit line. This is akin to the behavior observed with a linear resonator, but with a 
very different dependence on Pp (Schuster et al, 2005). First, an abrupt discontinuity 
in the ac-Stark shift at Qo/ric = 3.1 clearly indicates the sudden increase in the intra- 
cavity average photon number n as the resonator switches from L to H. In this region, 
two spectroscopic peaks are observed at a given pump power. Second, the linewidth 
narrows down at large n, in strong contrast with the linear resonator case where it 
increases linearly (Schuster et al., 2005 ). 

Our measurements of the qubit frequency shift Awgi versus Pp allow us to quan- 
titatively test the theoretical predictions made in section 1.4.1 For that purpose, 

P from Awoi using values for the Sf and 
and insetl. We then fit 



1.14 



we use Eq. (1.32a) to determine h = ja^ 
Kf calculated for our qubit parameters [see Fig. 
the resulting experimental n{Pp) curves with the n values calculated from the sole 
dynamics of the resonator (with the qubit in |0)) given by the square modulus of 
the solutions of Eq. ( 1.15[ ) using K and the measurement line attenuation A as the 
only fitting parameters. The agreement is excellent over the whole (J7,Pp) range for 
K/2tt = -625 ± 15 kHz and A = 110.8 ± 0.2 dB, which is consistent with the design 
value K/2n = —750 ± 250 kHz and with independent measurements of the line at- 
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Fig. 1.14 (a) Experimental (open symbols) average photon number n and corresponding 
fits (lines - see text) as a function of Pp for the same datasets Qo/Qc ~ 3.1 (circles) and 
Qq/Qc = 0.7 (triangles) as in Fig. Experimental points are obtained by converting the 
measured Auoi according to the ac-Stark shift model of the inset (solid line - see text ; the 
dashed line shows the linear ac-Stark shift approximation), (b) Qubit dephasing rate r2 mea- 
sured for the same datasets (open symbols) and calculated either by numerical integration of 



the system master equation (solid symbols) or using Eq. (1.32fl (solid lines). The horizontal 
dashed-dotted line indicates the intrinsic dephasing in zero cavity field. The measured para- 
metric power gain (dashed line) is also shown for comparison for flo/Qc ~ 0.7. (c) Complex 



cavity field amplitude ao (solid and dashed lines) calculated from Eq. (1.151 for increasing 
pump powers and for the same datasets. Both ao (squares) and ai (triangles) are also shown 
at six points labelled A-F. The separation between the two pointer states (segments) is to be 
compared to the uncertainty disk of a coherent state (open circles or disks shown at point C 
and at points A-F in panel b). 



tenuation A = 111 ± 2 dB. This demonstrates that measuring the ac-Stark shift of a 
qubit is a sensitive method for probing the field inside a nonlinear resonator and for 
characterizing its Kerr nonlinearity ( Castellanos-Beltran and Lehnert, 20071. 

Having obtained from the ac-Stark shift an accurate calibration of all the sample 
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parameters including the microwave power reaching the sample, we are able to test 
quantitatively our predictions regarding the qubit dephasing rate induced by pumping 
the nonlinear resonator. We show in Fig. 1.14 b) the measured dephasing rate r2 = ttw 



of the qubit, with w the full width at half maximum (FWHM) of a Lorentzian fit to 
the experimental data. In addition to the measurement-induced dephasing F^m, r2 
also includes a constant contribution F2 due to all other dephasing processes (mainly 
energy relaxation, flux noise and radiative broadening). For the sake of comparison, 
we also show the independently measured gain G of the KR operated as a parametric 
amplifier on a small additional signal frequency-shifted by 100 kHz from ujp/2'iT. Below 
and above He, ^2{Pp) peaks where G is maximum and near the bifurcation threshold 
respectively, and tends towards F^ at large n. 

We also show the predictions of Eq. ( 1.32f| calculated without any adjustable 
parameter. The qualitative features of our data (the peak of the dephasing rate near 
the maximum of the gain or the bifurcation threshold) are well reproduced. The direct 
correlation between dephasing rate and distance in phase space between the nonlinear 



resonator pointer states (distinguishability D) is graphically evidenced in Fig. 1.14 



Despite this qualitative agreement, quantitative discrepancies are observed when 
F2 is large. To understand their origin, we have calculated F2 by numerical integration 
of the full master equation Eq. (1.14) for the sample parameters and damping rates 
determined independently [see full symbols in Fig. 1.14 b)]. The agreement with the 
experimental data is excellent with no adjustable parameter over all parameters range, 
in contrast with experiments with a flux-qubit in which an unexplained shortening of 



the qubit relaxation time is observed when the KR is driven above bifurcation ( Picot 



et ai, 2008 


Serban et al., 2010 


1.1 


and the predictions of Eq. ( 


1.32f 



2010). This implies that the discrepancies between our data 



samples parameters, but only to the breakdown of the hypothesis I? <C 1 needed to 
derive this analytical formula, which is indeed only satisfied by our data when r2 
is below ~ 10 MHz. We thus confirm experimentally the theoretical predictions of 
section 1.4.1 in the parameter range where the parametric gain is negligible (G ~ 1) 
so that Eq. (1.32a) is valid and the pointer states are coherent states, the backaction of 
the driven KNR on the qubit is close to the quantum limit. Further work is needed to 
study the large gain regime, and in particular to establish whether the qubit readout 
by a JBA discussed in section [l.3| is quantum-limited or not. 



1.5 Conclusion 

In this chapter we have investigated experimentally and theoretically the coupling of 
a transmon qubit to a resonator with a Kerr nonlinearity, which constitutes one of the 
simplest extensions to the well-studied Jaynes-Cummings model. We have shown how 
effects arising from the resonator nonlinearity (bistability, parametric amplification, 
squeezing) qualitatively modify the dispersive qubit-resonator coupling. This allowed 
us in particular to obtain a qubit readout which has both a high fidelity and the poten- 
tial of being QND. The quantum backaction exerted by the intraresonator field onto 
the qubit was also studied. We stress that despite the much more complex system 
dynamics, we obtain qualitative agreement with analytical calculations and quanti- 
tative agreement with numerical simulations, indicating that the sample parameters 
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are controlled as well as in linear circuit QED. This opens the way to more elaborate 
experiments in which for instance the qubit could be used as a probe of the quantum 
fluctuations of the nonlinear resonator. In that way detailed tests of the theory of 
fluctuating nonlinear resonator could be performed. 
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